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^ , Abstract 

. The notion of chaotic behavior is examined for particle production in branch- 
J> I ing processes. Two types of branching are considered: non-Abelian gauge inter- 
^ ' action and an Abelian cascade model. Properties of the production processes 
. are investigated by Monte Carlo stimulation. The "temporal" behavior is stud- 
CD i ied by following the fluctuations in the multiplicities of each generation as the 
' branching evolves. The "spatial" behavior is described in terms of the fluctua- 
Q"^ ■ tions of the normalized factorial moments from event to event. The information 
dimension and a new entropy index are determined. When all the measures 
Q^l are taken together, they collectively give a strong suggestion that the QCD 
Q^' branching process is chaotic, while the Abelian cascade model is not. 
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. 1 Introduction 

Recently we have reported on the results of an investigation in the possible signatures 
of chaos in branching processes [Q. The aim is to determine whether the nonlinear, 
non-Abelian dynamics of the quantum Yang-Mills field possesses chaotic behavior. 
Since in such a dynamical system the number of degrees of freedom increases with 
time evolution, when the notion of time is not even well defined in the branching 
process, new measures of trajectory, distance, entropy, etc. must be introduced. We 
have found that the perturbative QCD branching shows signs indicative of chaos, 
whereas a model lacking the characteristics of QCD does not. In this paper we give 
the details of our study. 

In the case of classical non-Abelian dynamics special simplifying conditions that 
reduce the equations of motion to manageable size have been considered and chaotic 
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solutions have been found [0, |] . A more complete investigation of the gauge equations 
has to be done on the lattice, and it has been shown that the classical non-Abelian 
gauge theory generally exhibits deterministic chaos, whose Lyapunov exponents can 
be numerically computed [4-6]. 

The extension of the above classical problem to the quantum theory of Yang- Mills 
field is extremely difficult. The current state of knowledge about quantum chaos 
hovers round semi-classical problems in which classical trajectories are generalized to 
waves 1^ ||. That is totally inadequate for treating quantum fields, the proliferation 
of which in a collision process involves issues that are untouched in solving wave equa- 
tions. The first step toward formulating a feasible program to attack the problem is 
to dissociate the complications of nonperturbative QCD from the quantum dynamics 
of non-Abelian fields. In the perturbative domain where the QCD coupling is small, 
the nature of nonlinear, non-Abelian dynamics is fully present so that the signature 
of chaotic behavior in gauge theory should nevertheless show up in hard processes. 
That is an important step of simplification that renders the problem manageable. 
More specifically, one can narrow the scope and focus on branching processes with 
QCD splitting functions. The issue then becomes the search for measures that can 
reveal chaoticity in branching processes. 

In such a search it is necessary to keep in mind the special features of QCD. 
There are many processes that involve branching in real time, such as cell division in 
biological systems. The emission of photons by an accelerated charge can be regarded 
as a branching process with sequential ordering. However, for gluons reproducing 
gluons in the pure gauge theory, the time variable plays no role in the description 
of a state with n gluons whose momenta are precisely specified. In multiparticle 
production one works with momenta (or their variants such as rapidity), since they 
are what can be measured in a collision process. Then in the momentum space 
the concept of trajectory for a system of increasing number of gluons becomes ill 
defined. Without trajectories the notion of distance between trajectories is untenable, 
and the conventional method of examining chaos in classical nonlinear problems is 
inapplicable. It will therefore be necessary for us to look for other quantities that can 
describe the difference between states, which evolve from what correspond to nearby 
initial conditions in the classical problem. 

The major difficulty with the study of chaotic behavior of nonlinear gauge dy- 
namics is that the gauge fields are not directly measurable and that the branching 
processes cannot be tracked experimentally. In collisions at high energies, where the 
question of chaos in gauge dynamics arises, only the particles in the final states can be 
measured. Thus to verify any theoretical predictions about chaotic behavior, the loss 
of information at the end of the branching process must be quantified and presented 
in a form suitable for experimental determination. To that end we shall introduce an 
entropy index /Xg, which is measurable, and describes the degree of fiuctuation of the 
final particles from event to event. 

Before going into details, it may be helpful to state the general idea underlying 
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the work. In classical dynamics if the coordinates qa{t) and momenta Pa{t) specify 
the system, then the trajectory in phase space is well defined in the familiar way. In 
classical field theory the fields, Ea{x,t) and Ba{x,t) say, form a field configuration 
defined over all space x at any given time, and the change of the configuration as time 
evolves specifies a trajectory in a generalized sense. For quantized fields one works 
in the Fock space so that the number of quanta becomes a variable that specifies an 
important aspect of the state of the system. The distance between two trajectories in 
that case must therefore involve, among other quantities, the difference between the 
numbers of quanta in the system. Thus in particle production at high energy where 
the multiplicity n of produced particles is an essential observable, the fluctuation of n 
from event to event must be viewed as the consequence of swings of trajectories that 
have almost the same initial condition. With that connection in mind it is natural 
to generalize the conventional treatment of nonlinear dynamics to an approach that 
places emphasis on tracking the multiplicity fluctuations from event to event in high 
energy collisions. After describing the QCD dynamics of branching and a simple x 
model in Sec. 2, we treat the temporal behavior of the branching processes in Sec. 3. 
Then in Sees. 4 and 5 we consider the spatial properties of the final state. Concluding 
remarks are given at the end. 



2 Gluon Branching and the % Model 

To carry out this investigation it is necessary to use computer simulation to generate 
events of particle production through branching. Only then is it possible to study 
the nature of fluctuations of the final-state particles. For the simulation we shall 
use algorithms based on two opposite types of dynamics. One is, of course, the 
QCD dynamics. For simplicity we shall focus on only the pure gauge theory without 
quarks. At the end of the evolution we shall identify the partons as particles to avoid 
the complication of hadronization, which is inessential to the question of whether 
the gluon branching dynamics is chaotic. The other is a cascade model, to be called 
the X model, which has none of the features of the QCD; in particular, it does not 
have infrared and collinear divergences. It is studied in order to provide a contrast 
to the gauge theory so that our measures for chaoticity can be tested on these two 
contrasting branching dynamics. 

In both branching processes the initial parton has virtuality Q^, and successive 
branchings continue until the virtualities of all partons are < Ql- In pure-gauge QCD 
the splitting function at each vertex of branching is (for g gg) 

P(z)=6\^ + ^ + zil-z)] , (1) 
. z 1 — z J 

where z is the momentum fraction of the daughter parton in the frame where the 
mother parton's momentum is 1. In the x model we keep only the last term of ([^) 

P{z) = 6z{l- z) (2) 
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so that it has no divergences at z = and 1. The presence of those divergences is, of 
course, the source of comphcation for QCD that must be treated carefully. 

We follow Odorico's procedure to develop the algorithm for Monte Carlo sim- 
ulation of parton shower in QCD. Because of the soft gluon and coUinear divergences 
many partons are emitted at small z and small angles but are not resolvable. The 
probability that a parton at = can survive without emitting a resolvable parton 
until = Ql is the Sudakov form factor |P, |T0| 

' dt' as{t') /■i--o(i') 



A(Q,Q^) =exp 



dzPiz) 



(3) 



iQl V 2tx Jzoit') 

where the limits of the z integration define what is meant by resolvable. The first 
step of the simulation is to set zq = Qq/Q"^ and calculate from A(Q^, Ql) whether a 
parton starting from evolves to without branching. If not, then the value of 
at which the branching occurs is determined by solving 



where is a random number between and 1. With that g^, the value of z is then 
generated in accordance to the distribution P{z) given in ([^) keeping only what is in 
the range Q^/q^ < z < 1 — Q\lq^. With that 2;, the daughter partons are assigned 
the maximum virtualities t\ = zq^ and t2 = {1 — z)q^ , from which further evolutions 
are carried out by repeating the above procedure starting with ti and ^2- When the 
virtuality of a descendant reaches g^ < then that branch of the tree terminates. 
We shall put Qq = 1 GeV, to be definite. In the simulation we shall set = 3, 
Nf = so that 

a,(g2)=4vr/lllog(gVA2) , (5) 

where will be set at 250 MeV. 

We have also tried to simulate the parton shower using the algorithm of Weber |TD 
and found that the result differs very little from that based on the above algorithm 
of Odorico for the type of measures we calculate. Since Odorico's method is far 
more efficient, we have chosen to use it throughout this work. 

In the X model, because of no infrared and collinear divergences, there is no 
automatic g^ evolution. However, in order to compare its properties with those of 
QCD, we introduce by hand dependence of the branching process by requiring 
that at each vertex the two daughter partons have virtualities zq'^ and (1 — z)g^, 
when the mother virtuality is g^, and z is generated by using @ for < 2; < 1. 
As with QCD, we require branching to continue successively until the virtualities of 
all partons become < Q^. Since g^ is not degraded along a parton line, there are 
far more particles produced in the x model than in QCD for the same Q^, but that 
is immaterial, since our measure of chaoticity will not depend on Q/Qo- Clearly, 
the dynamics of branching in the x model is very different from QCD. We use it to 
exemplify the Abelian dynamics that has no infrared and collinear divergences. 
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3 Temporal Behaviors of Branching 



Each branching process can be represented by a tree diagram, since recombination 
of partons is not considered. The vertices of the tree could be ordered vertically in 
accordance to the values of of the mother partons. Then there are many diagrams 
with the same topology that can describe different evolution processes leading to the 
same number of particles at the end, but belonging to different final states. If we 
defer the considerations of the momenta of the final particles until the next section, 
we can simplify the problem by focusing on only the topology of the tree diagram. 
In that way we study the fluctuation of the particle multiplicity at the expense of 
ignoring their momenta. What is gained is the possibility of defining a trajectory in 
the multiplicity space for a branching process. 

In emphasizing the topology of a diagram let us draw all diagrams in such a way 
that all partons of the same generation are placed at the same level regardless of 
their q^. An example is shown in Fig. 1. Where the partons are horizontally has no 
significance in that diag] am. Vertically, the branching points of the same generation 
are placed at the same level. All partons that reach the final state (represented by 
the dashed line) have < Q^. No information concerning > Qq is carried by the 
vertical position of the lines in the diagram. Let i denote the generation, staring with 
i = for the initial parton at q^ — Q^. Let 6j denote the number of branching points 
at the ith generation. In Fig. 1 we have 6j = 1, 2, 3, 3, 1 for i = 0, • • • , 4. Wc can define 
a vector b = (6o, ■ ■ ■) with as many components as there are nonvanishing 6j. Thus 
for the diagram in Fig. 1 we have b = (1,2,3,3, 1). If this description were to be 
applied to an electron radiating photons, then a bremsstrahlung diagram would have 
b = (1,1,1,---). On the other hand, for cell reproduction where every cell subdivides 
into two, then b = (1,2,4,8, • • •). The description can be further simplified, if we 
define 

Xi = log2 h (6) 

and the corresponding vector x = {xq, ■ ■ ■). Since the minimum Xi is and the max- 
imum Xi is i, the two extreme vectors are x = (0, 0, 0, • • •) which is like bremsstrahlung, 
and X = (0, 1, 2, 3, • • •) which is like cell-reproduction. They are shown in Fig. 2 by 
the thick lines, where Xi is plotted against i. All possible tree diagrams of branching 
processes are represented by a line in between the two boundary (thick) lines, as 
illustrated by the thin line in Fig. 2. Such a line specifies a trajectory. 

For finite the branching cannot go without end, so there is always a maximum 
i, which we denote by imax- The value of imax varies from event to event, even if 
is fixed for all events. Two paths with the same imax may differ for < i < i^ax , 
and the total numbers of partons produced would be different. In general, for any 
two trajectories x and x' one can define di = \xi — x'j^\ at every i where Xi and x'^ 
exist simultaneously. Obviously, di can be regarded as the distance between the two 
trajectories. Since two trajectories in x space can cross, d^ may vanish at nonzero i. 
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Since i is an integer, Xi has discrete values. Thus the trajectories cannot be dense in 
the continuum space of i and Xi, but for the discrete support of the trajectories we 
may regard the collection of all possible trajectories as being dense in the sense that 
a trajectory can pass through any allowed value of {i,Xi). 

To relate the description even more closely to the classical treatment of chaos, 
we now address the question of sensitivity to initial condition. In classical nonlinear 
dynamics one can specify the point in phase space where a trajectory begins and 
vary that point n in as small a neighborhood A^^ as one chooses. In non-Abelian 
gauge dynamics or in the x model, the initial condition is that the initial parton has 
virtuality = Q^. Being quantum mechanical, it is not necessary to vary in a 
small neighborhood of Q^. Quantum fluctuation is sufficient to guarantee that the 
final state of the branching process will vary from event to event, even if all events 
start out at precisely the same Q^. Thus the relationship to the classical consideration 
of initial condition is as follows. Consider J\f random points in the neighborhood of 
the initial point in phase space of a classical trajectory. Any such point can be mapped 
to the beginning of a particular branching process. A set of Af events in the Monte 
Carlo simulation all starting out at the same then correspond to J\f trajectories 
beginning in the neighborhood A^^ in the classical case. With this correspondence 
between the two problems it is then sensible to suggest that instead of studying the 
distance between two neighboring trajectories, one should consider all N" events in 
the Monte Carlo simulation and examine the mean deviation from the average of 
the parton multiplicities. Indeed, the study of fluctuations of observable quantities 
will be the main theme of our approach to analyzing chaotic behavior of branching 
processes. 

In Fig. 3 we show schematically several possible trajectories of the branching 
processes, all started out at the same Q^. Not only can trajectories that have the 
same imax be different, those with different imax are even more dissimilar. For the 
purpose of studying "temporal" behavior, which corresponds to the dependence on 
the generation i in the branching process, only trajectories with the same i^ax can 
be compared. Thus it is necessary to know the distribution P (imax) of imax for 
all events. Fig. 4 shows the result of 10^ simulations for each of the two type of 
branching processes discussed in the previous section. The initial virtualities are such 
that Q/Qo = 10^ for QCD and 15 for the x model. These values are chosen so that 
the average multiplicities are comparable, as we shall show below. Both P {imax) are 
Gaussian-like distributions, with the width being larger for QCD. This is the first 
indication of more fluctuation for QCD as compared to the x model. In choosing a 
narrow band of imax for further investigation, we select the shaded regions in Fig. 4 
that are situated at the maximum of P (imax)- 

As mentioned before, time is not a wcll-dcfincd variable in the branching process. 
It is deemphasized when we focus on the topololgy of the tree diagrams. The gener- 
ation label i carries only a rough notion of time, since it is well known, as in family 
genealogy, that a son can be younger than a grandson. However, there is one quantity 
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that increases monotonically with i; it is the parton muhiphcity Ui 



i-1 



i-1 



1 + 2"^^' 



(7) 



j=0 



j=0 



even though Xi may rise and fall with i. Define the average parton multiplicity (rij) 
at generation i by 



where is the value of rii at i for the eth event and the sum is over all Af events 
having the same imax chosen. This quantity (rii) may be taken to play the role of 
time in the branching process, although no linear dependence is implied. In Fig. 5 
we show the simulated results on (n^) as functions of i for various values of Q/Qo- In 
both types of branching processes (rii) increases monotonically with i, the rate being 
fastest at midrange of i. The values of Q/Qo are chosen such that {rii) are roughly 
in the same range for the two types. 

For the measure that describes the fiuctuation of the trajectories, and therefore 
of the multiplicities, we use the normalized variance 

{n^) 

which clearly gives a measure of the average "distance" between trajectories. It is 
always positive, as a distance function should be. In Fig. 6 is shown the results of 
simulation for the two cases. The sustained increase in the lower half of i range for 
QCD is a distinctive feature that is not shared by the x model. 

A more transparent way of presenting this result is to plot Vi vs (rii), as shown in 
Fig. 7. It corresponds to plotting the distance d{t) against t in the classical problem, 
except that here (rii) is not exactly t, but some representation of it. Now, the behavior 
for QCD appears universal, i.e., independent of Q/Qo over a wide range. In the log- 
log plot the behavior is approximately linear. The same is not true for the x model: 
Vi increases initially, then drops precipitously, as (rii) is increased. Furthermore, 
the maximum decreases with increasing Q/Qo- Thus the uncertainty in the parton 
multiplicity increases with (rii) in the former case, but saturates and then decreases 
in the latter case. If that uncertainty can be regarded as a measure of chaos, then the 
QCD dynamics is chaotic, while the x model is not. We can express the dependence 
of Vi on {ui) in the QCD case in the form of a power law 




e=l 



(8) 




(10) 
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we may compare it to the definition of Lyapunov exponent A in d{t) ~ e^*. If In(nj) can 
be interpreted to correspond to t, then k may be said to play the role of the Lyapunov 
exponent. However, nothing firm can be proven about these correspondences, so the 
above remarks should only be taken as a possible orientation in interpreting the 
implication of the result obtained. 

From what has so far been done, it is not possible for us to find a criterion on 
the magnitude of n to signify robust chaoticity. Furthermore, (rij) and Vi are aspects 
of the parton state at generation i before the completion of the branching processes. 
Thus they are not experimentally measurable. For particle production in high energy 
collisions, the temporal behavior of branching is primarily a theoretical problem. For 
experimental verification of any theoretical prediction, it is necessary to focus on the 
final state, which is the subject of the next two sections. 



4 Factorial Moments, Entropy, and Information 
Dimension 

In the final state the complete information about the partons, now identified as ob- 
servable particles, is registered by their momenta. If the branching trajectories never 
change, then there is certainty in where to find those particles in the momentum 
space, and the entropy (to be defined below) is zero, corresponding to no loss of in- 
formation. However, quantum fluctuation alone is sufficient to cause fluctuation in 
the final-state momenta, so entropy is not expected to vanish. In multifractal analysis 
of complex patterns, the information dimension is a compact way of summarizing the 
dependence of the entropy on resolution scale Since multifractal analysis of the 



multiplicity fluctuation in high-energy collisions has been proposed [Il2|, |13| and car- 
ried out O, |l5i , it is natural to examine the information dimension for the branching 



processes considered here 

If pj denotes the fraction of particles in an event that fall into the jth bin of size S, 
then in the limit of many such bins in the system under study, the entropy is deflned 

by 

5 = -Ep.lnp, , (12) 



and the information dimension is 16 



Di = - \imSnn5 . (13) 

The study of self-similar behavior of particle production has been greatly facilitated by 
the use of normalized factorial momenta Fg, suggested by Bialas and Peschanski ||17|| , 
who flrst showed that the statistical fluctuations are flltered out by those moments. 
However, since F„ is deflned for integer g > 2, it cannot be used to determine the 
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information dimension, which involves the derivative of Fq with repect to g at g = 1. 
For that reason the Gq moments for all real g, positive and negative, were introduced 



for analyzing the multifractal structure of particle production data [0, 13, 15|. Their 
drawback is that the statistical fluctuations must be taken out of Gq explicitly by 
hand. More recently, a method for continuing Fq to noninteger g, while maintaining 
its virtue of not being contaminated by statistical fluctuations, has been developed 
[|T^ . We shall therefore make use of it to determine Di. 

Before describing the details of the calculations, let us first make clear the space in 
which the self-similar behavior is examined. With Zi denoting the momentum fraction 
of an ith generation daughter parton at a branching vertex, where the mother parton 
has momentum 1, the momentum fraction of a final particle is then 



X 



\{^^ ■ (14) 



where the product is taken over all generations of a particular path in the branching 
tree, leading from the initial parton to the final particle under consideration. Since all 
Zi are known in a simulated event, the values of x for all particles can be calculated; 
they are all crowded in the x = region, since their sum is 1. This is especially true 
for QCD branching because of the soft gluon divergence in (|l]). If we expand the x = 
region by using the variable C = — log^g ^i S^t a broad Gaussian-like distribution, 
as shown is Fig. 8, with the peak at C = 3 corresponding to x = 10~^. The situation 
is not as bad for the x model, but still most particles are found in the small x region, 
as shown in Fig. 9. This highly uneven distribution of x is inappropriate for data 
analysis that involves the partition of the space into small bins. 

One can define another variable X, in terms of which the distribution is much 



smoother [|T9|. Let the inclusive distribution in x, averaged over many events, be 



p(x). Then the cumulative variable X is defined by 

X{x) = j p{x) dx' / 1 p(x') dx' , (15) 

where xi and X2 are two extreme points in the distribution p(x), between which X 
varies from to 1. In terms of X the inclusive distribution p{X) is constant. For 
QCD branching we shall take p(x') to be the C, distribution as in Fig. 8; for the x 
model p(x') will be just the x distribution as in Fig. 9. An illustration of the p{X) 
distribution for QCD is shown is Fig. 10 after simulating 5 x 10'^ events aX Q/Qo = 10^. 
Note the expanded scale along the vertical axis. 

Now, in X space we divide the interval < X < 1 into M bins of width 6 = 1/M. 
Let rij be the number of final particles in the jth bin in any given event. The factorial 
moment of gth order is 

M 

/,(M)=M-i5:n,.K-l)---(n,-g + l) . (16) 
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The normalized factorial moment after averaging over all events is 



(/.) / ihr 



(17) 



Since (/g) is the average of nj\/{nj — q)\ over all bins and over all events, we may 
write it as 



if,) = E 



n=q 



{n - q) 



;i8) 



where P„ is the multiplicity distribution in a bin. In general, P„ may be expressed 
as a convolution of the statistical {S) and the dynamical {D) contributions to the 
multiplicity fluctuations: 



S^D 



dt— D{t) 



(19) 



where a Poissonian distribution has been used for S. Putting (O) into (|T8|) yields 



(/.) 



dtf'Dit) 



(20) 



This is the standard moment of the dynamical distribution that has no statistical 
contamination If D{t) = 5{t — n), so that P„ is Poissonian, then (fq) = n'^ and 
Fq = 1. Any deviation of Fq from 1 is an indication of the presence of dynamical 
fluctuation. 

Phenomenologically, one does not have access to D{t). Using P„ as input from 
experiment or from simulation, one can determine (fq), but only for integer q. One 
may replace the factorials in (|1^) by gamma functions in order to continue ([T8|) to 
noninteger q But that procedure does not result in Pg = 1 for Poisson distribution 
for all q. In Ref. ||18[ a continuation procedure is developed that guarantees Fq = 1 
for all q when D{t) is trivial. That method will not be described here, but will be 
used to determine Fg for our branching processes. 

For notational brevity let us write the RHS of (^) as it'*) Then we have 



{it/{t)Dr)D 



{{pMY) 



D 



(21) 



where p is the fraction of event multiplicity in a bin, as in ([T2|) . Thus from (|l^) we 
get 



S = -M(plnp) 



D 



InM- — lnP„ 
dq 



(22) 



9=1 



If Fq has a power-law behavior near q = 1, i.e.. 



Fq oc M"^' 



(23) 
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it then follows from ( |I3| ) that in the large M limit 



d(fq 



dq 



(24) 

q=l 



Note that this relationship is derived under the assumption that ( |22D is meaning- 
ful, which demands that Fq is well defined in a range of q around q = 1 and that 
dFq/dq contains only dynamical information without spurious contribution arising 
from improper continuation procedure [Q. 

If dynamically the bin multiplicity is evenly distributed among all M bins so that 



D{t) = 6{t — to), then Fq = 1 for all q (including non-integer q) and ( ^2]) gives 
S = InM. On the other hand, if dynamically only a fraction a of the M bins are 
populated, the rest empty, i.e. 

D{t) = aS{t~ti) + (l- a)5{t) , (25) 

then (fq) = ati and Fq = a^"'^ for all q. Hence, (^) yields S = In(aM). If, in 
particular, only one bin is non-empty, i.e., a = 1/M, then S = 0. Thus (pSj) offers a 
simple way of seeing how the entropy increases from to InM, as the fraction of non- 
empty bins increases from 1/M to 1. If there is only one non-empty bin, the dynamics 
is like one with a classical trajectory that has no dispersion in where the particles are 
to be found in the momentum space. Without loss of information the entropy is zero. 
If there is maximum dispersion corresponding to all bins being equally likely to be 
populated, then the entropy is maximum. 

We now apply this method of analysis to our simulated branching processes. First, 
we show in Fig. 11 the multiplicity distributions in QCD for various bin sizes (5 = 
1/M) after 10^ events. It should be understood that Fig. 10 shows the inclusive X 
distribution after summing over all events. But event by event the fluctuation in 
X space is much greater. Upon dividing the interval < X < 1 into M bins and 
counting the multiplicity in each bin, one then obtains P„(M) after averaging over 
all the events. The degree of fluctuation around the decreasing (n), relative to that 
average, as M is increased, is what is measured by the normalized factorial moments 
Fq defined by (|T^ and (|T^). This is unlike our theoretical calculation based on ( PP]) 
and (^) because D{t) is not known explicitly in the simulation. The results are shown 
in Fig. 12 for both QCD and x model. Since Fq=i = 1 by definition, we see that Fq 
increases with q in QCD, but decreases with increasing q in the x model. It means 
that the distributions Pn{M) for QCD are wider than Poisson, but for the x model 
they are sub-Poissonian. Fig. 12 also shows that there is very little M dependence 
in either case. For non-integer q values a substantial amount of work is needed to 
continue Fq according to the method of [0. The result is shown in Fig. 13. It is 
clear that for g > 1 all Fq are greater in QCD than in the x model. It corresponds to 
more fluctuations in QCD than in the x model, a result that is consistent with the 
finding of the previous section. 
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It is also evident from Fig. 13 that dFq/dq\^^^ is essentially independent of M in 
both cases. Whereas S is still dependent on M according to (0), the M independence 
of dhiFql dq\^^^ requires, on account of p3|), that dipq/dq\^^^ = 0. This together with 
(H) entails 

Di = 1 (26) 

in both cases. Thus there is no interesting multifractal property. This particular 
approach to finding a useful signature of chaotic behavior ends in failure. 

The smoothing of the phase space by use of the X variable in place of the mo- 
mentum fraction variable x has undoubtedly contributed to eliminating some aspect 
of the fractal structure. But what has not been removed is the more interesting con- 
sequences of the dynamics. The triviality of Di states only that there is no nontrivial 
self-similarity in the fluctuations measured. The latter are quantified by Fg, which 
involves both the bin and event averages of the factorial product in (|16|). The result 
suggests that some of the fluctuations have been lost by the averaging process. In the 
next section we develop a scheme to recapture that which has been lost. 



5 Entropy Index 

We know that event by event there are large fluctuations in where the produced 
particles are in the X space. What we want to do now is to find a way to register the 
fluctuation (from event to event) of the fluctuation (in the X distribution). 
Let F^ denote the normalized factorial moment of the eth event: 

Fq = n/{ny , (27) 

where is as defined in (|16|) for the eth event. Clearly, F^ changes with e for fixed 
q and M. In Fig. 14 we show the distribution P{Fq) after 10^ simulated events of 
QCD branching, the label e being omitted. For clarity, only q = 2 and 3 are exhibited 
for a wide range of M values. At small M the distribution has a peak at Fq ^ 1. 
Indeed, if M = 1 and event multiplicity is large, then Fq should be very nearly 1 for 
all events. At large M a peak at = is developed, especially at high q, because 
the bin multiplicity for small 6 is small and may in many events be less than q for 
all bins, in which cases F^ = . This large fluctuation in P{Fq) is what we want to 
capture, and would be lost if Fq is averaged over all events. Clearly we should take 
the moments of Fq before performing the average. 

We may think of F^ as the horizontal moments in the X space, and then define 
the vertical moments in the event space as follows: 

{F^) = lr1:{F^Y ' (28) 

■'^ e=l 
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where M is the total number of events and p is a positive real number not restricted 
to integers. The normalized moments of moments are 

CUM) = {f^,{M)) I {F,{M)Y . (29) 

If Cp^q{M) has a power-law behavior in M, i.e., 

CUM) oc M^'^?') , (30) 

then we can define 

^J^q= ^i'q{p)\p=l ■ (31) 

which will be referred to as the entropy index. 

To see the connection between /Xg and entropy, let us define 

Pq = Fql E Fq (32) 

e=l 

and a (new) entropy in the event space 

Sq = -Y.P^,\nP^, . (33) 

e=l 

Furthermore, define the moments 

^p,. = E(nT , (34) 

e=l 

which are related to Sq by 

Sq = -■^\\iHUp=i ■ (35) 
On the other hand, Hp g can be related to Cp g by 

Cp,q=X^'-'Hp,q , (36) 

SO we have 

d 

dp 

In then follows from (|30|) and ( pT)) that 

Sq = ln{ArM~^'^) , (38) 



■\nCUp=i = ^^^f-Sq . (37) 
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apart from a possible additive term that is independent of M and M. 

This entropy, as defined in (|33[) , is clearly different from the entropy 5* defined in 
(p!2D. To emphasize that Sq is defined in the event space, we may call it eventropy. 
We can think of the event space as a one dimensional space with M sites. At each 
site we can register a number F^. If is the same at each site, then = l/A/", 
and Sg = InM. We should think of this as being highly disordered in the event space, 
since F^ is spread out uniformly over all space. The larger the number of events, 
the larger is the eventropy. This is very similar to the situation where S = InM, 
when the bin multiplicity is uniformly distributed in the X space that has M bins. A 
branching dynamic that results in the same F^ for every event does not fluctuate in the 
branching processes. It corresponds to nearby trajectories staying nearby throughout. 
In short, the dynamics is not chaotic. For Sg to be InA/", we see from (^8]) that fiq 
must vanish. Thus small fig corresponds to large eventropy, which in turn implies no 
chaotic behavior. 

On the other hand, if we consider the other extreme where all F^ are zero except 
one event e', then = 6ee' and 5*^ = 0. This is highly ordered in the event space, 
but the fluctuation of F^ from zero to a nonzero value is large. Generally speaking, if 
the distribution P{Fg) is broad, the fluctuation is large, initially nearby trajectories 
become widely separated in the final states of different events, and the dynamics is 
chaotic. In order for the eventropy to be small, fig must be large. Thus large entropy 
index implies chaotic behavior. 

In Fig. 15 the results of our simulation are shown for g = 2, 3, 4 and p = 0.5, 1.0, 
1.5, and 2.0. From the log-log plots one can identify approximate scaling behavior. 
We used the region M = 5 — 20 to determine the value of /ip, for which Cp^g for an 
incremental region of p around p = 1 has to be examined. We obtain 

^(QCD) ^ 0.0061, 0.054, 0.23 (g = 2,3,4) , 
^(x) = 0.0014, 0.010, 0.046 (g = 2,3,4) . 

This result is plotted in Fig. 16. Clearly, /if is much larger than It is 

suggestive that the QCD dynamics is chaotic, while the x model is not. However, 
at this point we have no quantitative criterion on how large fig must be in order to 
be chaotic. In classical nonlinear problems there is also no criterion on how long the 
time lapse must be in order for the distance between nearby trajectories to become 
sufficiently far apart to qualify for divergent behavior. Nevertheless, the positivity 
of the Lyapunov exponent is a simple condition, the counterpart of which for the 
entropy index is lacking. 

6 Conclusion 

We have explored various ways of measuring chaoticity of branching processes. Be- 
cause the problem is far more complicated than a classical nonlinear problem, the 
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search has not been a straightforward extension of the conventional ideas. Our ap- 
proach has been based on a simple premise: if a branching process is chaotic, then 
it should be hard to predict what the likely outcome of an event would be, given the 
knowledge of the result of a previous event. It means that there are large fluctuations 
from event to event. 

To be more concrete, we have focused first on the temporal development of branch- 
ing. The variance Vi of the multiplicities at a particular generation i can be regarded 
as a measure of the mean distance between "trajectories". We have found that Vi 
increases with (n^) in a power-law fashion for QCD branching, but not for the x 
model, an example of Abelian branching dynamics. That is the first indication that 
the QCD dynamics has the characteristics of being chaotic. We have not investigated 
the behavior of the higher moments, Vi being related to only the second moment. It 
is probably safe to predict that the higher moments would increase even more rapidly 
with {rii). 

For fluctuations in the flnal state where experimental observations are feasible, 
we have considered Fg, Di, Cp^g, and /ig as possible measures of chaoticity. It turns 
out that the information dimension is not effective, but the entropy index can be. 
What emerges is a realization that one should study the phase space and the event 
space simultaneously. More specifically, consider a two-dimensional lattice, where the 
horizontal axis has M sites, corresponding to the M bins in the X space, and where the 
vertical axis has Af sites, corresponding to the Af events in the event space. For every 
event there is a horizontal array of M numbers, indicating the particle multiplicities 
in the M bins. The moments of moments Cp^g summarize the distribution on this 
whole lattice without significant loss of information. The entropy index /ig extracts 
from this lattice of numbers a simple numerical quantity that can characterize the 
degree of fluctuations on the lattice. The larger Hg is, the more chaotic the system 
is. Unfortunately, we do not have at this point a criterion on the threshold value of 
^q, above which chaotic behavior will definitely occur. Further investigation on this 
aspect of the problem is undoubtedly warranted. 

Each of the features studied reveals some aspect of the chaotic behavior of QCD 
dynamics. Taken together collectively, they present a strong contrast from the results 
of the X model. One may wonder what the mechanisms are in the branching processes 
that can lead to such diverse outcomes. Our view is that in QCD branching there is 
a tension between two opposing forces that is absent in the x model. The coUinear 
divergence implies a strong preference for small angle emission of partons in the 
branching process. But if the emission angle is too small to be resolved, there is 
no branching. Since evolution without branching is suppressed by the Sudakov form 
factor, there is a competition between evolution by emitting a resolvable gluon and a 
preference for emitting a collinear unresolvable gluon. It is this tension that leads to 
large fiuctuations in the number of branching points and hence to chaotic behavior. 
The X model has no collinear divergence and no evolution, so there is no tension to 
cause large fluctuations. What we learn here may lay the foundation for flnding ways 
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to treat nonperturbative processes, for which chaotic behavior may play a crucial role 
in the evolution of a large system of quarks and gluons, such as those produced in a 
heavy-ion collision at high energy. 
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Figure Captions 

Fig. 1 A tree diagram for b = (1, 2, 3, 3, 1) 

Fig. 2 A trajectory of Xi 

Fig. 3 Several possible trajectories. 

Fig. 4 Distribution of i^ax for QCD and the x model. 

Fig. 5 Simulated results for the average multiplicities at various generations i for 
fixed values of Q/Qo indicated by the numbers. 

Fig. 6 Simulated results for the normalized variances for the same parameters as 
in Fig. 5. 

Fig. 7 Temporal behaviors depicted by Vi vs (n,). 

Fig. 8 Momentum distribution in terms oi — log x for QCD branching. 
Fig. 9 Momentum distribution in terms of x for the x model. 

Fig. 10 Momentum distribution in terms of the cumulative variable X for QCD 
branching. 

Fig. 11 Multiplicity distributions in bins of size S = 1/M for QCD branching. 
Fig. 12 vs M for g = 2 - 5. 
Fig. 13 Fq vs continuous q. 

Fig. 14 Distributions of calculated event by event. 

Fig. 15 Moments of moments Cp^q vs M for p = 0.5 — 2.0 and g = 2, 3, 4. 
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